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First-principles calculations are used in order to investigate phonon anomalies in non-magnetic 
and magnetic Heusler alloys. Phonon dispersions for several systems in their cubic L2i structure 
were obtained along the [110] direction. We consider compounds which exhibit phonon instabilities 
and compare them with their stable counterparts. The analysis of the electronic structure allows 
us to identify the characteristic features leading to structural instabilities. The phonon dispersions 
of the unstable compounds show that, while the acoustic modes tend to soften, the optical modes 
disperse in a way which is significantly different from that of the stable structures. The optical modes 
that appear to disperse at anomalously low frequencies are Raman active, which is considered an 
indication of a stronger polarizability of the unstable systems. We show that phonon instability 
of the TA2 mode in Heusler alloys is driven by interaction(repulsion) with the low energy optical 
vibrations. The optical modes show their unusual behavior due to covalent interactions which are 
additional bonding features incommensurate with the dominating metallicity in Heusler compounds. 



I. INTRODUCTION 

In this paper, we present a detailed discussion of the 
electronic, magnetic and vibrational properties of Heusler 
compounds. This subject is of particular interest in re- 
lation to the magnetic shape-memory (MSM) effect re- 
cently discovered in these materials Q . The Heusler com- 
pounds are known since 1903. However, only in the last 
decade have promising applications been presented for 
those systems which undergo a structural phase trans- 
formation and a magnetic phase transition near room 
temperature. The origin of the structural instability is 
related to the physics of martensitic transformations, of 
which the driving forces are not completely understood 
on a microscopic level 0. Therefore, a detailed knowl- 
edge of the difference between structurally stable and un- 
stable Heusler compounds on a microscopic scale would 
lead to a significant contribution to the theory of marten- 
sitic transformations in general and would be helpful 
when discussing further technical applications of MSM 
alloys. 

With respect to technological applications, Heusler al- 
loys arc important for two reasons. Those which undergo 
a martensitic transformation can be used to develop me- 
chanical devices based on the specific elastic properties 
of the martensitic structure, whereas the ferromagnet- 
ically ordered martensites can be used for MSM tech- 
nology 0. The MSM physics is based on the fact that 
the ferromagnetic martensitic structure can be deformed 
by applying an external magnetic field, since in MSM 
systems, martensitic domains are at the same time mag- 
netic domains. It takes less energy to redistribute the 
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FIG. 1: (a) The L2i Heusler structure with unit cell of 
space group Fm3m(05) of Ni2MnGa and (b) conventional 
tetragonal cell used in the electronic structure calculations. 
All compounds considered in this work have the same struc- 
ture and formula unit X2YZ. 



martensitic domains in case of strong uniaxial magnetic 
anisotropy compared to the energy required to rotate the 
local axis of magnetization. Corresponding domains that 
have a favorable orientation of magnetization grow at 
the expense of domains for which the magnetization di- 
rections are not parallel to the external magnetic field. 
This has been experimentally observed in Ni2MnGa 0], 
which today serves as a reference system for all investiga- 
tions related to MSM technology Q. It has been shown 
that in moderate magnetic fields of the order of 1 T, the 
structural deformations in Ni2MnGa can reach 10% 0. 
This is the key feature of MSM technology that allows 
the design of new kinds of micro-mechanical sensors and 
actuators. Therefore, a complete understanding of the 
interplay of martensitic micro-mechanics and magnetism 
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is of prime technological interest. 

In this paper, we present results of investigations of 
the structural stability and electronic properties of the 
L2i ordered non-magnetic (NM). ferromagnetic (FM) 
and fcrri-magnetic (FcrriM) Heusler compounds (the L2i 
structure is shown in Fig. From this compara- 

tive study, we expect to obtain those critical parame- 
ters which distinguish unstable and stable Heusler com- 
pounds. The present analysis is a continuation of previ- 
ous investigations of structural instabilities in such com- 
pounds, whereby the instability has been related to Fermi 
surface nesting in FM compounds 0j H @1- In addi- 
tion to the nesting picture, we present in this paper fur- 
ther unusual features of Heusler alloys which are related 
to characteristics of phonon dispersions and atomic dis- 
placements, electron and phonon densities of states, and 
the coexistence of the metallic and covalcnt bondings. 



we add one valence electron in the system by substitut- 
ing Ge to the place of Ga(Al,In), and obtain compound 
Ni2MnGe with e/a = 7.75. In Ni2MnGa and Ni2MnGe, 
we substitute Ni by Co in order to consider compounds 
CosMnGa (e/a = 7.00) and CoaMnGe (e/a = 7.25). A 
crossover from the Ni-based to the Co-based compounds 
is especially important due to their very different proper- 
ties. The Co-based compounds are known as half-metals 
[11 [13. Eventually, we obtain more or less systematic 
sampling of different compounds which, from prelimi- 
nary analysis, are known to be ferromagnetic. In order 
to discuss the role of magnetic order [^ along with fer- 
romagnets we consider one ferri-magnetic and one non- 
magnetic systems, FeaMnGa and NiaTiGa, respectively. 

III. METHOD 



II. HEUSLER ALLOYS SELECTED FOR THIS 
STUDY 

A set of eight Heusler alloys have been selected for 
this investigation. They are connected by a close rela- 
tionship to NiaMnGa which is a reference system for all 
studies related to the MSM properties and Heusler alloys 
in general, exhibiting most of the properties of interest. 
Moreover, we are quite familiar with this system from the 
previous theoretical investigations [3, 0| ■ At compo- 
sitions close to stoichiometric, this compound has the 
cubic L2i structure (see Fig. which is stable at high 
temperature. Under cooling this crystal undergoes a se- 
ries of magnetic and structural phase transitions 
Of central interest is the soft phonon mode shown by the 
L2i structure of NiaMnGa at low temperature The 
character of this soft mode is believed to be related to 
the modulated structures 3M, lOM and 14M, coupling 
to a uniform tetragonal distortion and the specific nest- 
ing topology of the Fermi surface of NiaMnGa [^IMITol|. 
Similar properties have been found in other Heusler com- 
pounds [13 ■ By making controlled changes to the chemi- 
cal composition of NiaMnGa wc plan to explore chemical 
trends in the electronic structure and structural energet- 
ics. 

It can be helpful to characterize different Heusler com- 
pounds by their average concentration of valence elec- 
trons, which is defined as an electron-to-atom ratio 
[13I lig . [l7l | . The parameter e/a is a measure of the elec- 
tronic filling and allows in a rough way to distinguish 
in the phase diagram of the ternary alloys the range 
of stability for the martensitically unstable alloys. In 
case of NiaMnGa, Ni contributes ten, Mn seven, and Ga 
three valence electrons yielding e/a = 7.5. We consider 
those compounds which have the same number of va- 
lence electrons as NiaMnGa and those whose e/a are 
slightly different. Thus, by substituting one atom of 
Ga in NiaMnGa by atoms of Al or In, we obtain com- 
pounds NiaMnAl and NiaMnIn with e/a = 7.5. Then, 



A. First principles calculations 

The Vienna Ab-initio Simulation Package (VASP) 
[20I I2H was used to perform the first-principles calcula- 
tions. The projector-augmented wave formalism (PAW) 
implemented in this package [21I l2^ leads to very ac- 
curate results comparable to other all-electron meth- 
ods. The electronic exchange and correlation are treated 
within density functional theory by using the generalized 
gradient approximation {2^ . The expansion of the elec- 
tronic wave-functions in plane waves was done using the 
"High Precision" option in VASP, which corresponds to 
the kinetic energy cutoff as high as 337.3 eV or more, de- 
pending on the system. For the total-energy calculations 
in the conventional tetragonal unit cell shown in Fig^b), 
we used a high-density mesh of 12x12x10 special k- 
points for Brillouin-zone integrations. For the larger su- 
percell used in the phonon calculations, described below, 
we used a 10x2x8 mesh. 



B. Phonon calculations 

Phonon dispersions along [110] were calculated from 
first principles for each system in the L2i structure at 
the theoretical lattice constant given in TableQ] We used 
the direct force constant method [3[13- The supercell 
was a 1 X 5x 1 periodic supercell based on the conventional 
tetragonal cell with lattice parameters at = bt = a/\/2, 
Ct = c, shown in Fig^b). This is an orthorhombic super- 
cell with the long axis along [110], containing ten consec- 
utive (110) atomic planes along the [110] direction. We 
considered in turn each Cartesian displacement of each 
crystallographically independent atom by 0.03 A and 
computed the induced forces acting on all other atoms 
within the supercell using the Hellmann-Feynman theo- 
rem. From these results, we constructed the force con- 
stant matrix within the harmonic approximation. The 
phonon frequencies and corresponding eigenvectors were 
obtained by diagonalizing the corresponding dynamical 
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matrix. The 1x5x1 supercell geometry allows the direct 
calculation for five q- vectors along the [110] direction. 
These five vectors satisfy equation 

exp (27riqL • L) = 1, 

where L denotes indices of the lattice constants in the 
supercell (in our case from to 4, giving points — 
0.0, 0.25, 0.5, 0.75, 1.0, for the normalized wave vector 
[(^, C,0], which spans our Brillouin zone from its center 
to the boundary). The interpolation of the dispersion 
curves is performed by setting the force constants with 
range beyond five atomic planes (half the length of our 
supercell) to zero; in the present case the decay of force 
constants with distance is found to be quite rapid and 
thus this is an excellent approximation. 



IV. RESULTS 

A. L2i structure and magnetic ordering 

For the stoichiometric Heusler compounds (X2YZ) 
considered in this work, the calculated lattice constants, 
magnetic moments per formula unit and the nature of 
magnetic ordering in the L2i structure are given in Table 
^ All atoms occupy high symmetry positions within the 
Fni3m (^^225) symmetry, where WyckofF positions are 4a 
(0,0,0) for Y atoms; 4b (1/2, 1/2, 1/2) for Z atoms; 8c 
(1/4, 1/4, 1/4) for X atoms. The computed lattice con- 
stants largely agree with the experimental values, where 
available, as summarized in Ref. l26l l27i The magnetic 
moment of Fe2MnGa is the result of local magnetic mo- 
ments of Mn and Fe being aligned anti-ferrimagnetically. 
Further discussions of the magnetic pro perties of the 
Heusler alloys can be found in Ref. |23, 129 . 



B. Calculated phonon dispersions 

The calculated phonon dispersions of the eight selected 
Heusler compounds in their theoretical cubic L2i struc- 
tures are shown in Fig. [2| Some of the phonon disper- 
sions, shown in Fig. [3 (Ni2Mn(Ga,Al)), have already 
been discussed in recent publications |l5l l30l | . Compar- 
ison of the phonon dispersions of Ni2MnGa with exist- 
ing experimental data and with phonon dispersion cal- 
culations by other groups employing the linear response 
method show that the direct method used here yields 
fairly accurate results comparable to the former method 

Ultil. 

The [110] dispersions of the L2i structure have one set 
of nondegenerate acoustic branches (L, TAi and TA2) 
and three optical branches that can be easily recognized 
in Fig. 12 For the five compounds Ni2Mn(Ga,Al,In,Ge) 
and Ni2TiGa, the TA2 branch is unstable for some range 
of (. In addition, in Ni2MnGe the TA2 mode has a nega- 
tive slope at F, indicating a pure elastic instability. The 



instability of the L2i structure in NM Ni2TiGa shows 
that magnetic order is not a necessary condition for the 
phonon softening to occur. 

For some of the Heusler compounds the optical modes 
are split into three well separated groups which are triply 
degenerate at F, but become more mixed with increasing 
wave vector. In each of these three groups one finds dom- 
inating vibrations of a particular kind of atom. Mostly, 
this is due to the differences in the atomic masses and 
it can nicely be established by computing the partial 
phonon density of states for the [110] branch, shown for 
two selected compounds in Fig. 13 

Comparison of these two plots shows that the position 
of the Ni-peak in case of unstable compound Ni2MnGa is 
different from the position of the Co-peak in stable sys- 
tem Co2MnGc. This was unexpected, because atoms of 
Ni and Co, in both compounds, occupy the same sites in 
the structure, Wyckoff positions 8c(l/4, 1/4, 1/4). Their 
masses are close to each other. Thus, the sequence of the 
optical modes should be similar. However, while in case 
of stable compound Co2MnGe the sequence of the opti- 
cal modes is normal, for the case of Ni2MnGa we observe 
an inversion of the optical modes. 



V. DISCUSSIONS 
A. Inversion of the optical modes 

In the following we will discuss the impact of the in- 
version of the optical modes shown in Fig. |2| in relation 
to the traditionally accepted point of view of the nature 
of incommensurate instabilities in Heusler alloys. As a 
matter of fact the structural instability of the L2i phase 
is associated with Fermi surface nesting which also causes 
instability of the TA2 mode in some systems leading to an 
incommensurate product phase, which has a modulated 
structure. The relation between soft mode behavior and 
specific nesting feature of the Fermi surface has been re- 
cently emphasized and is not a debate any more 0, S 0] ■ 

However, a systematic review of the phonon disper- 
sions reveals a characteristic behavior of the optical T2g 
mode at F (see Fig. ^J. In all Ni-based systems we find 
that this mode appears at much lower energies as com- 
pared to the Co-based stable compounds. The role of this 
anomaly might be as important as the instability of the 
TA2 mode, or in other words, these two features have to 
be considered together in order to explain the structural 
instabilities in Heusler alloys. 

We drew attention to the fact that phonon modes of 
the same symmetry must repel each other. This is what 
must happen in case of the TA2 and the T2g modes. At 
finite wave vectors only one optical mode which has [1-10] 
polarization will repel with the the TA2. The main point 
is that there are two instabilities in the system. One of 
them is that the (110) planes in the Heusler structure 
can slide in the [1-10] direction, this is the TA2 mode. 
But at the same time the (111) planes of Ni slide against 
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TABLE I: Computed lattice parameters, magnetic moments per unit cell and types of magnetic order for the eight selected 
Heusler compounds in the L2i structure. 'Instability' of the cubic structure means here that an unstable mode appears in the 
calculated phonon dispersion. The valence-electron-to-atom ratio e/a for the NM compound is shown in brackets in order to 
avoid direct comparison with the magnetic systems in terms of this parameter. 



System 


aL2, (A) 


Mtotal (Mb) 


L2i 


Magn. order 


e/a 


Ni2MnGa 


5.8067 


4.35 


unstable 


FM 


7.50 


Ni2MnAl 


5.7000 


4.20 


unstable 


FM 


7.50 


Ni2MnIn 


6.0624 


4.22 


unstable 


FM 


7.50 


Ni2MnGe 


5.8039 


4.10 


unstable 


FM 


7.75 


Co2MnGa 


5.7100 


5.09 


stable 


FM 


7.00 


Co2MnGe 


5.7285 


4.99 


stable 


FM 


7.25 


Ni2TiGa 


5.8895 


0.00 


unstable 


NM 


(6.75) 


Fe2MnGa 


5.6882 


2.15 


stable 


FerriM 


6.50 



each other also along the same [1-10] direction, and this 
is the mode, shown in Fig. ^ These vibrations are 
destructive for each other and have to repulse. Thus, the 
TA2 mode is unstable because it is pushed down by the 
corresponding from the symmetry optical T2g mode of 
Ni. The Fermi surface nesting determines at which wave 
vector the TA2 mode is most sensitive to the influence of 
the optical mode. 



B. Hybridization features 

In order to understand this observation we have to re- 
call that Heusler alloys are very unusual metallic systems 
due to the presence of the so called p-elements (Ga,Al,Ge 
etc), which form partially filled bands close to Yjp. These 
bands allow the d electrons of transition metals hybridize 
with the p electrons of the p element j^^. In case of 
Ni2MnGa, the Ga atoms form energetically favorable hy- 
brid states with Ni. These states give a peak in the spin- 
down electronic density of states right at the Fermi level 
(see Fig. 13)). The presence of this peak around Ej? is 
related to the anomalous vibrational properties of the 
optical mode T2g- 

In contrast to this, Co2MnGe does not show a peak in 
the DOS close to Y,p. The same kind of hybridization 
probably occurs in this case, but the p band in this case 
sits at lower energies so that when d electrons of Go enter 
this band we observe what is called half-metallic behav- 
ior, i.e. there are no spin-down states at the Fermi level 
as can be seen in Fig. |S1 

This hybridization features affect binding mechanisms 
in Heusler alloys to a large extent. These intermetallics 
and a competition of covalent interaction and magnetic 
ordering was already discussed by Kiibler et al. in [28l |. 
In our case, the covalent bond stems from the Ap elec- 
trons of Ga which can couple to the 3d electrons of Ni. 
The magnetism is mainly governed by the magnetic Jahn- 
Teller effect. This effect splits the giant peak of Mn 3d 



states, which is present in the non-magnetic DOS at the 
Fermi level 0- The stabihty of the FM Heusler sys- 
tems depends on a balance between magnetic ordering 
of the Mn atoms and covalent bonding of Ga and Ni 
atoms. However, we have shown, taking Ni2TiGa as an 
example, that non-magnetic Heusler compounds exhibit 
the same kind of martensitic instability as the FM ones 
and the feature in DOS is also present in this case (see 
Fig. EJ. Therefore, the magnetic order cannot be con- 
sidered to be at the origin of the structural stability. In- 
stead, more general features of the band occupations and 
the hybridization have to be considered. It is important 
that while in metallic systems the covalent interactions 
are expected to be weak, the present calculations of the 
phonon dispersions show that some Heusler systems can 
exhibit strong anomalies due to the covalent bonding. 

It is important that the low lying optical modes of Ni 
are the only Raman-active modes. We know from the 
mechanism of Raman scattering that this process must 
involve the polarizability of the material. This is cer- 
tainly related to the covalent interaction between the Ga 
Ap and Ni 3d states. The distortion of the electron cloud 
will help us further to consider one more aspect of this 
issue because the tetrahedral environment of the atoms 
in the Heusler alloys leads to a natural distortion of the 
charge around X atoms (in X2YZ; X=Ni,Co,Fe). 

Apparently, it is difficult to obtain experimental evi- 
dence of the optical modes inversion since the metallic 
character of our systems will wash-out any signals of the 
Raman activity. However, Raman-active Ni vibrations 
have so far been observed in an amorphous metal |33l |. 
A Raman scattering experiment on Ni2MnGa would al- 
low to follow up the vibrational behavior of Ni from the 
parent phase through the martensitic structure with de- 
creasing temperature. We suggest that the anomalous 
sequence of the optical phonons in the structurally un- 
stable Heusler systems is a general feature and can be 
applied generally to other intermetallics. 

Finally, since we have involved the issue of optical 
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FIG. 2: Phonon dispersion curves of (a) FM Ni2MnGa, (b) FM Ni2MnAl, (c) FM Ni2MnGe, (d) FM Ni2MnIn (e) FM Co2MnGa, 
(f) FM Co2MnGe, (g) NM Ni2TiGa and (h) FerriM Fe2MnGa in the L2i structure. Here, the reduced wave vector coordinate ( 
spans the fee Brillouin zone from F to X. Imaginary frequencies of the unstable modes are shown in the real negative frequency 
range. The frequency of the optical modes T29 at F is marked by a black; note it appears in lower values as compared to the 
stable systems. 



phonons and covalent interactions in our discussions, it 
is worth thinking here about possible similarity with fer- 
roelectrics, in which the structural instability can also 
be discussed on the basis of hybridization effects, in this 
case of oxygen p and transition metal states. For the 
unusual behavior of phonon dispersion in ferroelectrics 
see, for example, |33 |. 



C. Trends with e/a 

The parameter e/a is of rather qualitative nature and 
has to be discussed in relation to other details of the elec- 
tronic structure discussed above. Here, we present some 
observations which can be useful in the development of a 
conceptual macroscopic view on the stability of Heusler 
structures. 

The character of softening of the acoustic TA2 mode of 
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FIG. 3: Site projected vibrational density of states showing 
the distribution of amphtudes of the atoms over the frequency 
range. This information shows which atom dominates in the 
vibrations of a given frequency range. In particular, the 
anomalous behavior of optical vibrations of Ni lying at fre- 
quencies below those of the the heavier Ga atom in unstable 
Ni2MnGa is clearly visible (see also Fig.l^J. In the case of the 
stable system Co2MnGe, the sequence of the optical modes is 
regular, i.e. increasing as the mass of the atom decreases. 



Ni2MnGa, Ni2MnAl and Ni2MnIn is similar. For these 
compounds e/a is the same, 7.5. If we compare the 
phonon dispersions of these systems to those having one 
additional valence electron, i.e. e/a = 7.75, as in the 
case of Ni2MnGe, we find that the onset of instability at 
C w 0.6 is more pronounced. We do not focus here on 
the additional instability in the dispersion of Ni2MnGe 
at r, it has a different origin. In contrast to this unstable 
behavior of the Ni-based systems, the phonon dispersions 
of ferromagnetic Co-based stoichiometric compounds do 
not show any trace of a structural instability (Fig.|2Jl. In 
this case, the number of valence electrons is smaller than 
7.5. Thus, we see that Co- and Ni- based systems show 
different behavior with respect to the shear instability in- 
volving the TA2 phonon mode. Although, it is not clear 




FIG. 4: (Color online) The F-point T29 optical mode which 
vibrates at anomalously low energy. Note that this is the only 
phonon mode with T2g symmetry in this structure and it is 
Raman active. Moreover, this mode is pure Ni-vibrations, of 
Co, i.e. the atoms sitting in the 8c Wyckoff positions. Other 
atoms are not involved in this mode. Other optical modes 
involve vibrations of all atoms together and have symmetry 
Tiu at F, being IR-active. 
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FIG. 5: (Color online) Site projected electronic density of 
states of Ni2MnGa and Co2MnGe in the L2i structure. There 
is a difference in the density of spin-down electronic states of 
both compounds close to the Fermi level. In Ni2MnX (X = 
Ga, Ge, Al, In) compounds, the peak in the DOS marked by 
the arrow arises from p-states of the X atom. The double 
structure of this peak is related to the hybridization splitting 
cause by the interaction of 4p| and Ni-3d| states [3^ . 
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FIG. 6: (Color online) Site projected electronic density of 
states of Ni2TiGa. The double peak structure at the Fermi 
level is caused by the hybridization of Ga 4p and Ni 3d states 
showing that non-magnetic Heusler systems can also exhibit 
the same kind of characteristic peak structure as the FM sys- 
tems. 



what is that key difference between Co and Ni, which 
makes Heusler structures stable and unstable. Intrinsic 
properties of the Ni and Co atoms would require different 
kind of approach not used in this paper. It is left for fu- 
ture discussions in order to concentrate on other aspects 
of our work. 

We come back to the value of e/a in order to check 
if there is a trend which suggests a critical e/a value to 
determine whether the L2i phase is stable or unstable. 
Maybe it is possible regardless what kind of atom (Ni 
of Co) we have in the structure, which would be very 
useful for analizing alloying effects in the Heusler alloys 
and their stability. For this purpose we consider nearest 
neighbor force constants between the atoms sitting in the 
8c Wyckoff positions i.e. Ni or Co in our case. Figure 
13 presents results for only ferromagnetic systems. From 
this picture we can see that both for Ni- and Co- based 
systems the force constants follow similar trends suggest- 
ing lost of their stability at some critical values for e/a. 
Here we have two possibilities. We extrapolated lines 
corresponding Co- and Ni- based systems and show that 
thay have different behavior. Such an approach shows 
that replacing Co by Ni would shift the linear depcn- 
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FIG. 7: (Color online) Nearest neighbor Ni-Ni and Co-Co 
force constants (X atoms in X2YZ Heusler compounds). The 
force constants strongly depend on the valence electron num- 
ber per atom. Possibly, there is a critical value of e/a allowing 
the separation of structurally unstable Heusler systems with 
negative force constants from those of stable ones with posi- 
tive force constants, but this has to be verified by additional 
calculations. Note that these Ni-Ni negative force constants 
do not lead directly to the unstable acoustic mode TA2, but 
it leads to the inversion of the optical modes. 



dence down and other way around if we repkace Ni by 
Co. Contrary to that, we can assume that the kind of 
atom is not important, but the e/a dependence of the 
force constants is not linear and a crossover from stable 
to unstable configurations occurs at around e/a = 7.4. 
However, our data is not sufficient to draw any final con- 
clusion for these two scenarios. This needs to be verified 
in future studies. 

We would like to note that when looking at three sys- 
tems, Ni2MnAl, Ni2MnGa and Ni2MnIn with the same 
e/a ~ 7.5 we find that the differences in their force con- 
stants are small. This shows that the size and mass of Al, 
Ga or In do not play much role, but the similar valence 
electron filling seems to be more important. 

We did not include Ni2TiGa and Fe2MnGa in the FigEI 
keeping in mind that the type of magnetic order should 
lead to significant distinctions in how the electrons fill 
appropriate energy levels with possible relation to the 
critical ratio of e/a. With respect to the role of the 
magnetic order, it is helpful to keep in mind previous 
results related to the interplay of magnetic order and 
martensitic instability has previously been discussed for 
Ni2MnGa It has been shown that the spin-split 

band structure gives rise to a prominent peak in the the 
wave vector dependent magnetic susceptibility close to 
q = ^(i,i,0), which corresponds to the experimen- 
tally observed wave vector for the prcmartensitic insta- 
bility 0. This enhanced susceptibility has its origin in 
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Fermi surface nesting in the different spin channels. The 
fact that FM Co-based systems are stable, although ten- 
dencies of Fermi surface nesting exist, reflects an obvious 
situation when the characteristic topology of the Fermi 
surface must strongly depend on the e/a ratio. 



D. Tetrahedral coordination of the covalent 
bonding 

In order to further elucidate the role of the covalent 
interaction, we have analyzed the distribution of the elec- 
tronic charges in our ferromagnetic Heusler alloys. The 
comparison of the spin-up and spin-down densities of 
states has shown a clear difference in the charge distribu- 
tion. The issue of the electronic hybridization becomes 
clear if we analyze the situation in Ni2MnGa by using a 
simple model. Figure |S1 shows schematically the charge 
distribution (spin-down only) taken from ah initio calcu- 
lations for Ni2MnGa in the (110) plane of the tetragonal 
unit cell shown in Fig. ^p. We observe an electrostatic- 
like repulsion between the Mn and Ni atoms. The reason 
for this effect is the large magnetic moments of Mn. Spin- 
down 3d electrons of Ni do not find symmetry-allowed d- 
states of Mn to hybridize with, unless they flip their spin. 
Instead, the id electrons of Ni hybridize with the Ap elec- 
trons of Ga. It looks like the charge of Ni is pushed to- 
wards the neighboring Ga atoms, which are tetrahedrally 
coordinated around it. Thus, the tetrahedral coordina- 
tion which we flnd in the Heusler alloys is a natural basis 
for the id electrons of Ni to hybridize with the Ap elec- 
trons of Ga. 



VI. CONCLUSIONS 

Our total energy density functional theory investiga- 
tions of the structural instabilities in non-magnetic and 
magnetic Heusler compounds allow us to draw following 
conclusions. 

We have found that in all Heusler compounds which 
show acoustic unstable mode TA2, the optical vibrations 
exhibit unusual inversion of their modes. Analysis of the 
force constants revealed that negative force constants are 
present in all unstable systems, but these force constants 
are not directly related to the acoustic phonon instabil- 
ity. Instead, the inversion of the optical modes appears 
to be the driving force for the acoustic anomaly. The 
Raman-active optical modes lower their energy due to 
the negative force constants, whereby the acoustic mode 
TA2 becomes unstable due the repulsive interaction with 
the anomalous optical mode. This acoustic-optical inter- 
action might be of more general interest then just for the 
Heusler compounds. Its role in the martensitic transfor- 
mations has to be investigated in more detail, including 
those in binary alloys. 

A comparative analysis of different systems has shown 
that unstable Heusler compounds show characteristic fea- 




FIG. 8: (Color online) Schematic representation of the defor- 
mation of the outer 3d electron cloud around the Ni atoms 
in the unstable Heusler alloys from present ah initio calcu- 
lations. There are two effects which contribute to the defor- 
mation. One is due to covalent bonding of Ga 4p and Ni 3d 
states, the other (stronger) effect is due electrostatic repul- 
sion of Mn and Ni atoms resulting from their nearly spherical 
charge distribution. The deformation of the 3d orbitals of Ni 
is of tetrahedral nature and strengthens the covalent Ni-Ga 
bond. 



tures in the electronic DOS. The peak in the electronic 
states appears at the Fermi level, being responsible for 
the nesting topology of the Fermi surface of the unstable 
Heusler systems and the covalent bonding features which 
affect the optical vibrational modes of the Ni-based sys- 
tems. 

We obtained that phonon instabilities observed in the 
Heusler compounds can be related to the number of va- 
lence electrons in the systems. However, our consider- 
ation of only Co- and Ni-based systems does not allow 
us to discuss how general this observation can be. A 
crossover from the stable Co- to unstable Ni-based sys- 
tems on the basis of the e/a value can be suggested, but 
additional investigations are required. Microscopic in- 
trinsic properties of Co and Ni atoms may stand for a 
substantial part of our observations. 

Study of additional Heusler compounds (see Ref. l27j) 
will serve to further establish physical trends to the prob- 
lem of martensitic phase transitions. 
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